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Abstract 

In a two-way relay network, two terminals exchange information over a shared wireless half-duplex 
channel with the help of a relay. Due to its fundamental and practical importance, there has been an 
increasing interest in this channel. However, surprisingly, there has been little work that characterizes the 
fundamental tradeoff between the communication reliability and transmission rate across all signal-to- 
noise ratios. In this paper, we consider amplify-and-forward (AF) two-way relaying due to its simplicity. 
We first derive the random coding error exponent for the link in each direction. From the exponent 
expression, the capacity and cutoff rate for each link are also deduced. We then put forth the notion 
of the bottleneck error exponent, which is the worst exponent decay between the two links, to give 
us insight into the fundamental tradeoff between the rate pair and information-exchange reliability of 
the two terminals. As applications of the error exponent analysis, we present two optimal resource 
allocations to maximize the bottleneck error exponent: i) the optimal rate allocation under a sum-rate 
constraint and its closed-form quasi-optimal solution that requires only knowledge of the capacity and 
cutoff rate of each link; and ii) the optimal power allocation under a total power constraint, which 
is formulated as a quasi-convex optimization problem. Numerical results verify our analysis and the 
effectiveness of the optimal rate and power allocations in maximizing the bottleneck error exponent. 
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I. Introduction 

The two-way communication channel was first introduced by Shannon, showing how to 
efficiently design message structures to enable simultaneous bidirectional communication at the 
highest possible data rates [1]. Recently, this model has regained significant interest by introduc- 
ing an additional relay to support the exchange of information between the two communicating 
terminals. The attractive feature of this two-way relay model is that it can compensate the spectral 
inefficiency of one-way relaying under a half-duplex constraint [2]-[7]. With one-way relaying, 
we should use four phases to exchange information between two terminals via a half-duplex 
relay, i.e., it takes two phases to send information from one terminal to the other terminal and 
two phases for the reverse direction (see Fig. [T]). However, exploiting the knowledge of terminals' 
own transmitted signals and the broadcast nature of the wireless medium, we can improve the 
spectral efficiency by using only two phases to exchange information in the two-way relay 
channel (TWRC) [2]. 

Due to the aforementioned fundamental and practical importance of the TWRC, much work 
has investigated the sum rate and the achievable rate region of the TWRC with different re- 
laying protocols [2]-[7]. The half-duplex amplify-and-forward (AF) and decode-and-forward 
(DF) TWRCs have been studied in [2] where it was shown that both protocols with two-way 
relaying can redeem a significant portion of the half-duplex loss. In [3], the achievable rates 
for AF, DF, joint-DF, and denoise-and-forward relaying have been analyzed and the condition 
for maximization of the two-way rate are investigated for each relaying scheme. The broadcast 
capacity region in terms of the maximal probability of error has been derived in [5] for the 
DF TWRC. A new achievable rate region for the TWRC has been found in [6] for partial DF 
relaying, which is a superposition of both DF and compress-and-forward relaying. Bit error 
probability at each terminal has also been analyzed for a memoryless additive white Gaussian 
noise (AWGN) TWRC [7]. However, there has been few work that characterizes the fundamental 
tradeoff between the communication reliability and transmission rate in the TWRC across all 
signal-to-noise ratio (SNR) regimes. 

In this paper, we consider half-duplex AF two-way relaying due to its simplicity in practical 
implementation. To characterize the fundamental tradeoff between the communication reliability 
and rate, we first derive Gallager's random coding error exponent (RCEE) — the classical lower 
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bound to Shannon's reliability function (see, e.g., [8]-[ll] and references therein) — for the link 
of each direction in the AF TWRcQ Instead of considering only the achievable rate or error 
probability as a performance measure, the RCEE results can reveal the inherent tradeoff between 
these measures to unveil the effectiveness of two-way relaying in redeeming a significant portion 
of the half-duplex loss in the information exchange. From the exponent expression, the capacity 
and cutoff rate for each link in the TWRC are further deduced. We then introduce the bottleneck 
error exponent, which is defined by the worst exponent decay between the links of two directions, 
to capture the tradeoff between the rate pair of both links and the reliability of information 
exchange at such a rate pair. Using this notion, we can appertain a bottleneck exponent value to 
each rate pair and characterize the bottleneck exponent plane from the set of all possible rate 
pairs besides the achievable rate region. This enables us to design a two-way relay network with 
reliable information exchange. 

For applications of the error exponent analysis for the TWRC, we present two optimal resource 
(rate and power) allocations, the main results of which can be summarized as follows. 

• We show that the optimal rate allocation to maximize the bottleneck error exponent under a 
sum-rate constraint is a rate pair such that the RCEE values of both links become identical 
at the respective rates. This optimal rate pair can be determined by a closed-form solution 
for sum rates less than a certain constant — called the decisive sum rate — depending only on 
the cutoff and critical rates of each link. Furthermore, the optimal solution requires only the 
knowledge of each cutoff rate. At sum rates larger than the decisive point, we can allocate 
a rate pair quasi- optimally in closed form, requiring only knowledge of the capacity and 
cutoff rate of each link. 

• We determine the optimal power allocation that maximizes the bottleneck error exponent 
under a total power constraint of the two terminals. In the presence of perfect global channel 
state information (CSI), we show that this power allocation problem can be formulated as a 
quasi-convex optimization problem, where the optimal solution can be efficiently determined 
via a sequence of convex feasibility problems in the form of second-order cone programs 
(SOCPs). 

The rest of this paper is organized as follows. In Section II, we describe the system model. In 

' In the following, we shall use simply the term "TWRC" to denote the AF TWRC. 
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Section III, we present the results of the error exponent analysis for the TWRC. The optimization 
framework for two-way relay networks is developed for the rate and power allocations to 
maximize the bottleneck error exponent in Section IV. We provide some numerical results in 
Section V and finally conclude the paper in Section VI. 

Notation: Throughout the paper, we shall use the following notation. Boldface upper- and 
lower-case letters denote matrices and column vectors, respectively. The superscript (■)^ denotes 
the transpose. We use M, M+, and to denote the set of real numbers, nonnegative real 
numbers, and positive real numbers, respectively. A circularly symmetric complex Gaussian dis- 
tribution with mean fi and variance cr^ is denoted by CM (fi, a^) and the exponential distribution 
with a hazard rate A is denoted by £ (X). 



We consider the TWRC as illustrated in Fig. [H where a half-duplex relay node R bidirec- 
tionally communicates between two terminals TkeT={i,2} with AF relaying. In the first multiple 
access phase, the terminals TkeT transmit their information to the relay and the received signal 
at the relay is given by 



CM (0, is the channel coefficient from to the relay, and 2r ~ CM (0, A^o) is the complex 
AWGnI] Note that l/i^f ~ ^ (V^^a,-)- 

At the relay, the received signal is scaled and broadcasted to both terminals in the second 
broadcast phase, while satisfying its power constraint pr. Then, the received signal at the terminal 
Tfcer is given by 



II. System Model 



Vr = hiXi + h2X2 + 2;r 



(1) 



where Xk(zr is the transmitted signal from the terminal Tker with E||xfcp} = Pk , 



yT,k = hkX-R + 2;T,fc 



(2) 



where x-^ = Gyn is the transmitted signal from the relay with E{ Ixr^} = pr, ZT,k ~ CM (0, Nq) 
is the AWGN, and G is the relaying gain given by 




(3) 



^ We assume the channel reciprocity for /ifc as in [2]. 
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We impose a total transmit power constraint P such that pi+p2 < P- As in [2], we further assume 
that the terminal Tker knows its own transmitted signal and has perfect CSI to remove self- 
interference prior to decoding. For notational convenience, we shall refer to the communication 
link Ti ^ R ^ T2 as the link Li and T2 ^ R ^ Ti as the link L2, respectively. With the 
self-interference cancellation, the effective SNR of the link Lf^^j- is given by 

Pfcttfc + {pr + Pm/Pk) aia2/ak + 1 

where a^er = l/ifelV^o- 
III. Error Exponent Analysis 
A. Mathematical Preliminaries 

The reliability function or error exponent for a channel of the capacity C is the best exponent 
decay with the codeword length N in the average probability of error that one can achieve at a 
rate R<C [8]-[ll]: 

E (R) = lim sup In P^* {R, N) (5) 

where P°p* (R^N) is the average block error probability for the optimal block code of length 
and rate Ra As a classical lower bound on the reliability function, the RCEE or Gallager's 
exponent is given by [8] 

Ej- (R) = max max {Eo{p, Q )—pR] (6) 
Q 0<p<l 



with 



Eo(p,Q) = -ln 



dy } (7) 



where Q (x) is the input distribution and p{y\x) is the transition probability. Unfortunately, 
the double maximization in ^ is generally very difficult since the inner integral is raised to 
a fractional exponent when p E (0, 1) and the lack of knowledge about the optimal input 
distribution Q (x). For analytical tractability, the Gaussian input distribution Q{x) is often used, 
which is optimal if the rate R approaches the channel capacity [8]-[ll]. 



' Throughout the paper, we shall use a rate measured in units of nats per second per Hz (nats/s/Hz). 
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B. Two-way Relay Channels 

The TWRC consists of two communication links and the achievable rate can be characterized 
by the sum rate of two parallel relay channels under perfect self-interference cancellation [2]. 
As such, we need to first consider the RCEE for each link and subsequently introduce a notion 
of the bottleneck error exponent for the TWRC to effectively capture the tradeoff between the 
individual rates and the reliability. Using the Gaussian input distribution, we obtain the following 
proposition for the RCEE of each link in the TWRC. 

Proposition 1: With the Gaussian input distribution, the RCEE for the link Lfcgr of the TWRC 
with AF relaying is given by 



where 



E,^k (R) = max {Eo,fc (p) - 2pR} (8) 

pe[o,i] 



Eo,fc (p) = - InE^eff <( ( 1 + ^ ] }. (9) 



1 + P 



Proof: It follows immediately from the results of [10] along with the self -interference 
cancellation at the terminal TkeT- □ 
Remark 1: It is difficult to obtain a closed-form solution for Q in Proposition [U due to an 
analytically intractable form of the effective SNR 7^1^ given in (|4]). In what follows, to alleviate 
such difficulty and render ^ more amenable to further analysis, we use the upper bound 7^'^ 
on the effective SNR 7^^^^ by ignoring the term 1 in the denominator: 

^ub _ PfcPR«ia2 ^^Q^ 

PkOik + {pr + Pm/Pk) aia2/ak 
which corresponds to the ideal/hypothetical AF relaying [12]-[14]. 

Remark 2: The factor 2 of pR in ([8]) is due to the use of two phases for the exchange 

of information in the TWRC. In contrast, with one-way relaying, the information exchange 

occurs over four phases and hence, this factor should be 4, leading the RCEE for each link to 

Er,k (R) = maXpe[o,i] {Eo,k (p) - 4pi?}I] 

^ In the one-way relay channel (OWRC), if the total relaying power for information exchange is again constrained to pa, 
then the ideal/hypothetical AF relaying yields the upper bound on the effective SNR for the link L^gr as 



Pk (pr/2) aia2 



PkOik + (pr/2) a-ia2/ak 
which is slightly different from l IlOl l but makes no deviation in the analysis. 
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Theorem 1: With the Gaussian input distribution, the RCEE for the link L^gr of the TWRC 
with ideal/hypothetical AF relaying is given by 



E,,k{R) = max \Eo,k{p)-2pR} 
pe[o,i] I J 



(11) 



with 

Eo,k (p) 



-In 



4:XkPk 



y/nT (p) (v^+ y^fc) 



_ tt1,1,1,1,2 
4-"l,(l:l),0,(l:2) 



Vk 



(l+p)(v^+v/Mfc) 



(2,1) 

(l-p,l); (1/2,1) 
(0,1) ;(0, 1),(0,1) 



2 (Afe + pfc) 



v^vrr (p) (v^+ ^/Pk) 



4-'-'l,(l:l),0,(l:2) 



Vk 



(l+p)(v^+v^)' 
(V^fc+V^) 



(2,1) 

(l-p,l); (1/2,1) 

(0,1); (1,1), (-1,1) 
for < p < 1 (12) 



tK,N,N',M,M' 



and Eo^kip) = for p = 0, where F (■) is Euler's gamma function, p ^^.j^^ 

generalized Fox i/-function [15, eq (2.2.1)], and 

PR 



is the 



Vk 



PR+PlP2/Pk 

No 
Pk^k 



{PlP2/Pk + Pr) ^1^2 

Proof: See Appendix |Al □ 
Remark 3: The maximum of the exponent E^ ^ (R) over p occurs at i? = | dEo^k (p) / dp 
and hence, the slope of the exponent-rate curve at a rate R is equal to — 2popt. The maximizing 
Popt lies in [0, 1] if 



P=Popt 



'cr,k 



dEp^k (p) 
dp 



< R< 



p=i 



dEp^k (p) 
dp 



iCk) 



(13) 



p=0 



where -Rcr.fc and (C^) are the critical rate and the (ergodic) capacity for the link Lker, respectively. 
For R < RcT,k, we have popt = 1, yielding the slope of the exponent-rate curve is equal to 
—2 and E^^k (R) = Ep^k (1) — 2i?. Furthermore, the cutoff rate for the link L^er is given by 
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i?o,fc = k (1) /2. This quantity is equal to the value of R at which the exponent becomes zero 
by setting p = 1. While the capacity determines the maximum achievable rate, the cutoff rate 
determines the maximum practical transmission rate for possible sequential decoding strategies 
and indicates both the values of the zero-rate exponent and the rate regime in which the error 
probability can be made arbitrarily small by increasing the codeword length. 

Corollary 1: The ergodic capacity for the link Lfcgr of the TWRC with ideal/hypothetical AF 
relaying is given by 



_ 7^-1,2,1,1,2 
4-"l,(2:l),0,(2:2) 



(Afc + /Xfc) \/Xkf^k 0-1,2,1,1,2 

T?-" 1,(2:1), 0,(2:2) 



Vk 



(\/^+v^) 



Vk 
W>^kf^k 



(2,1) 

(1,1), (1,1); (1/2,1) 

(1,1), (0, 1); (0,1), (0,1) 
(2,1) 

(1,1), (1,1); (1/2,1) 



(1,1), (0,1); (1,1), (-1,1) 



(14) 



Proof: See Appendix |Bl □ 
Corollary 2: The cutoff rate for the link L^^r of the TWRC with ideal/hypothetical AF 
relaying is given by 



0,fc 



■iln 

2 



_ 17-1,1,1,1,2 
4-'-'l,(l:l),0,(l:2) 



2 (Afc + /ifc) yXkfJ'k 7^1,1,1,1,2 

v4-"l, (1:1), 0,(1:2) 



Vk 



2(v^+v^)' 
(v^+V^fc) 



Vk 

2(VA^+VMfc) 

W>^kf^k 
(v^+v^) 



(2,1) 

(0,1); (1/2,1) 

(0,1) ;(0, 1),(0,1) 
(2,1) 

(0,1); (1/2,1) 
(0,1); (1,1), (-1,1) 



(15) 

□ 



Proof: It follows immediately from (|T2l) by setting p = I. 
Remark 4: It is insufficient to characterize the information exchange in the TWRC by only 
investigating the RCEE for each link individually, as it just reflects the tradeoff between the 
communication rate and reliability for the information transmission in one direction. Therefore, 
we introduce a notion of the bottleneck exponent for the TWRC to capture the tradeoff between 
the rate pair of both links and the reliability of information exchange at such a rate pair, enabling 
us to optimize the resource allocation in the TWRC. 
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Definition 1 (Bottleneck Error Probability): For a TWRC with the terminal T^gr transmitting 
a code (e^^*, A^) of rate Rk, the bottleneck error probability is defined as 

P* = maxPj*^) (16) 

where Pc^^ is the error probability of the link h^. 

Note that Definition [T] can be applicable for a general TWRC, regardless of relaying protocols. 
From the random coding bound 

p{k) < g-A^^,,fc(Hfc) (^^Y^ 

the bottleneck error probability of the TWRC is bounded by 

P,* < maxe-^^^''=(^*\ (18) 
~ keT 

Using (fTSi) . we define the bottleneck error exponent of the TWRC as follows. 

Definition 2 (Bottleneck Error Exponent): For a TWRC with the terminal T^er transmitting 
a code (e^^'^jA^) of rate Rk, the bottleneck error exponent at the information-exchange rate 
pair (Ri, R2) is defined as 

E,*(Pi,P2) =minE,,fc(Pfc). (19) 



Remark 5: Using the RCEE of the link Lker in Theorem [H we can readily obtain the 
bottleneck error exponent E* (Pi, P2). From Definition [2l we can see that the bottleneck error 
exponent captures the behavior of the worst exponent decay between the two links in the TWRC 
and reflects the reliability of the information exchange at a rate pair (Pi, P2). When the worst 
link is good enough, it means that the other link must also be good. As a result, using (fT9l) as 
an information-exchange reliability measure, we can design a two-way relay network such that 
both links can communicate reliably. Besides the achievable rate region, we can also characterize 
the bottleneck exponent plane from the set of all possible rate pairs. This plane could provide 
us with further understanding of the tradeoff between the rate pair (Pi, P2) and the bottleneck 
error exponent (i.e., information-exchange reliability). 

IV. Optimal Resource Allocation 

A. Optimal Rate Allocation 

In the following, we present the optimal rate allocation that maximizes the bottleneck error 
exponent P*(Pi,P2) under a sum-rate constraint in the reliable information-exchange region 
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^ = -R2) : < -Ri < (Ci), < R2 < (C2)}. Mathematically, this rate allocation problem 
can be formulated as follows: 



^1 



max E*{Ri,R2) 

Ri,R2 



S.t. 



(20) 



Ri + R2 = R 

0<Ri< (Ci), < i?2 < {C2) 
which can be solved by the following theorem. 

Theorem 2: Let ^ and ^ be the curve and straight line in ^ such that 

^ = {(i?i,i?2) e ^ : ^,,1 = E,,2 (i?2) 

if = i?2) G : i?i + i?2 = R} . 

Then, the optimal solution (i?i,i?2)opt of the rate allocation problem Vi for the sum rate R > 
|i?o,i — -Ro,2| is the intersection point of the rate-pair curve ^ and straight line In particular, 
we have 

/ R + -Ro,i — Ro,2 R — Rqa + Ro,2 



i?2 



'opt 



(R,0) 
(0, R) 



for |i?o,i -^0,2! < R < Rd 
for R < i?o,i — Ro,2, Ro,i > Ro,2 
for R < Rq 2 — Ro,i, Ro.i < Ro,2 



where R^ is the decisive sum rate given by 

Rd = '^i'^ {2-Rcr,l — -Ro,l + -Ro,2, 2-Rcr,2 + -Ro,l " -^0,2} • 

Proof: See Appendix O 
Remark 6: For the sum rate R > Rj, we can determine the quasi-optimal rate pair as 



(21) 

(22) 

□ 



/ R + -Ro,l — Rq,2 R — Rq,! + -Ro,2 



R2 



(23) 



2 > 2 , for R^ < R < R^ 

'opt - WR - (C2), (C2)) for R > R^, (Ci) > {C2) 

^ ((Ci), R - (Ci)) for R > R^, (Ci) < (C2) 

where R^ = min {2(Ci) — i?o,i + -Ro,27 2(C2) + -Ro.i — -Ro,2}S Therefore, with knowing the ca- 
pacity and cutoff rate of each link in the TWRC, we can determine the optimal rate pair 

^ The numerical example in Section |V] will show that this quasi-optimal rate pair well approximates the optimal one for the 
sum rate R > RJ. 
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(i?i,i?2)opt '^hat maximizes the reliability of information exchange at the sum rate R — exactly 
for R < using (|2TI) . and approximately for R > R^ using (l23l) . 



B. Optimal Power Allocation 

In this subsection, we present the optimal power allocation that maximizes the bottleneck error 
exponent E* (Ri, R2) at a rate pair R2). In the presence of perfect global CSI, for fixed p 
and (_Ri,i?2), we are maximizing the instantaneous bottleneck exponent over p = [pi p2]^ for 
each fading state, i.e., before averaging with respect to fading. Mathematically, we can formulate 
this optimization problem as follows: 



max El^'ip,p,R,,R2) 
p 



s.t. 



Pl+P2<P 
Pi > 0, P2 > 

where the subscript "int" denotes an instantaneous value and 



(24) 



El^'ip,p,R,,R2) = mm El^lip,p,R,) 



KlijP.P.Rk 



In 



1 + 



PkPRaia2 



(25) 



- 2pRk. (26) 



1 + ppkQk + {pr + pm/Pk) aia2/ak + 1. 
With the optimizing p^^^ obtained by solving the problem (|24l) . we can find the bottleneck error 
exponent with optimal power allocation as follows: 



E: R2) = E,,,„, ( max E'l'^' {p^^„ p, R^, R2) ] . 



(27) 



Since pi and p2 are positive, we can define ipk&T — JPk and -0 = \}\)\ ip2]^ without loss of 



optimality. With this change of variables, we can transform the optimization problem in (1241) 
into a quasi-concave programl^ 

Theorem 3: For fixed p and rate pair (_Ri, R2), the function E™* (p, p, R2) is quasi-concave 
and the program V2 is quasi-concave. 

Proof: See Appendix |Dl □ 



* Let 5 be a convex subset of R^. A function / : 5 R is said to be quasi-convex if and only if its lower-level sets 
L {f,a) = {x £ S : f (x) < a} are convex sets for every a G R. Similarly, / is said to be quasi-concave if and only if its 
upper-level sets U {f,a) = {x £ S : f (x) > a} are convex sets for every a £ R. 
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It is well known that we can solve quasi-convex optimization problems efficiently through a 
sequence of convex feasibility problems using the bisection method [16]0 We formalize it in 
the following corollary. 

Corollary 3: The program V2 can be solved numerically using the bisection method: 
0. Initialize tmin and tmax, where tmin and tmax define a range of relevant values of 
-^r* (P' P' -^1' -^2), and set the tolerance e G 

Solve the convex feasibility program Psocp (t) in (l28l) by fixing t = (tmax + ^min) /2. 



1. 
2. 
3. 
4. 



t else set 



t. 



If S (t) = 0, then set t, 

Stop if the gap (tmax — i^^min) is less than the tolerance e. Go to Step 1 otherwise. 
Output ■0opt obtained from solving Psocp (t) in Step 1. 
where the convex feasibility program can be written in SOCP form as [17] 

^socp (t) : find ip 

s.t. tp eS{t) 

with the set S (t) given by 

tjj'^ei/y/vl ] [ ip^e2/y/v^ 

Sit) = {^e 



(28) 



Alp 

a/1 + PR«2 



0, 



a/1 +PRai 



^yc 0, 







(29) 



where A and tifegr are defined by (|52l ) and (|53l) in Appendix |Dl respectively. 

Proof: It follows directly from the proof of Theorem [3] and [17] that we can represent the 
convex constraints in the set S (t) in terms of SOC constraints. □ 
Remark 7: It is important that we initialize an interval that contains the optimal solution. In 
our case, we can always let tmin correspond to the uniform power allocation and we only need 
to choose tmax appropriately. 

V. Numerical Results 

In this section, we provide some numerical results to illustrate our analysis. In all examples, 
we choose = 0.5, = 2, pr = P, and define SNR = P/Nq. Without power allocation, we 
further consider equal power allocation between two terminals T^er, namely, pi = P2 = P/2@ 

^ Note that the program P2 is always feasible as long as P > 0. 

^ For one-way relaying, the RCEE, capacity, and cutoff rate are symmetric and equal for two links in the case of equal power 
allocation. 
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A. Random Coding Error Exponent 

To ascertain the effectiveness of two-way relaying in terms of the error exponent, Fig. [21 shows 
the RCEE for the link Lfc^r of the TWRC and OWRC with ideal/hypothetical AF relaying at 
SNR = 20 dB. To calculate the RCEE, we use Theorem [T] for two-way relaying, whereas we 
modify Theorem [U for one-way relaying in such a manner as described in Remark [21 We can 
see from the figure that the link L2 of the TWRC shows better exponent behavior than the link 
Li at every rate R due to the fact that VL2 > ^i- In the regime below the critical rate, the 
exponent of the TWRC decreases with the rate twice as slow as in the OWRC and hence, we 
require to increase the codeword length slowly with two-way relaying to achieve the same level 
of reliable information exchange as the rate increases. This is due to the spectral efficiency of 
two-way relaying that requires only half the time duration of one-way relaying to exchange the 
information. 

B. Capacity and Cutoff Rate 

Figs. [Hand m demonstrate the effectiveness of two-way relaying on the achievable rates, where 
the capacity (or achievable sum rate) and cutoff rate versus SNR are depicted for the link L^g-r 
of the TWRC and OWRC with ideal/hypothetical AF relaying, respectively. We can see from the 
figures that the slopes of the capacity, achievable sum rate, and cutoff rate curves at high SNR 
are twice as large in the TWRC as in the OWRC due to again the fact that two-way relaying 
for the information exchange can reduce the spectral efficiency loss of half-duplex signaling by 
half in the TWRC. Hence, as can be seen in Fig. [3l the high-SNR slope of the capacity for the 
link Lfcgr of the TWRC is identical to that of the achievable sum rate in the OWRC. 

C. Bottleneck Error Exponent 

To demonstrate the tradeoff between the rate pair and information-exchange reliability in the 
TWRC, Fig. [51 shows the bottleneck error exponent E* {Ri, R2) versus Ri for the TWRC with 
ideal/hypothetical AF relaying at SNR = 20 dB when R2 = 0, 0.2, 0.5, 0.7, and 1.1. For fixed 
R2, the bottleneck exponent (Ri, R2) as a function of -Ri behaves identically to Er,i{Ri) 
at Ri > Rf"" = mm{Ri G M+ : E,^i{Ri) < ^^,2(^2)}, whereas E^{Ri,R2) is limited to 
Er,i {Rf"") for all Ri < Rf"". In this example, the values of Rf"" are equal to 0, 0.16, 0.36, 
0.54, 0.92 for R2 = 0, 0.2, 0.5, 0.7, and 1.1, respectively. As can be seen, the bad link in terms 
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of the exponent is a bottleneck that limits reliable information exchange, and the bottleneck 
exponent at large Ri or becomes small, indicating the achievable reliability of information 
exchange would be low at such rate pairs. 

D. Optimal Resource Allocation 

We now give application examples of the error exponent analysis for the resource allocation 
in the TWRC. 

1) Optimal Rate Allocation: Fig. [6] shows the optimal rate pair (i?i,i?2)opt that maximizes 
the bottleneck error exponent E*{Ri,R2) under a sum rate constraint for the TWRC with 
ideal/hypothetical AF relaying at SNR = 20 dB. The quasi-optimal rate pairs are also plotted for 
the sum rate R > R^. The optimal and quasi-optimal rate pairs are determined using Theorem |2l 
and (|23l) . respectively. The optimal rate pairs (-Ri,-R2)opt '^he sum rates R = 0.148, 0.4, 0.8, 
1.2, 1.6, 2.0, and 2.4 are (0, 0.148), (0.126, 0.274), (0.326, 0.474), (0.520, 0.680), (0.710, 0.890), 
(0.910, 1.090), and (1.103, 1.297), attaining the maximum E* (i?i, R2) equal to 1.37, 1.12, 0.73, 
0.37, 0.15, 0.04, and 1.8 x 10""^, respectively. For R > R^ = 0.83, the quasi-optimal rate pairs at 
the sum rates R = 1.2, 1.6, 2.0, and 2.431 are (0.526,0.674), (0.726,0.874), (0.926, 1.074), and 
(1.118, 1.282), attaining E* {Ri, R2) equal to 0.366, 0.1449, 0.0316, and 0, respectively. We can 
see that the quasi-optimal rate pairs quite well approximate the optimal (-Ri, -R2)opt > ^^d 

and achieve the bottleneck exponents very close to the maximum achievable E* (Ri, R2) at such 
sum rates. In the figure, the region ^ can be divided by the optimal rate curve into two rate-pair 
subregions in which each RCEE Ej.^^ (Rk) is dominant for the bottleneck exponent E* (i?i, R2): 
for example, Ej-^i (Ri) is dominant in the light gray subregion, i.e., E* {Ri, R2) = E^^i (Ri)- 

The effectiveness of the optimal/quasi-optimal rate allocation in maximizing the bottleneck 
error exponent can be further ascertained by referring Fig. Ul where the bottleneck error exponent 
E* (i?i, R2) versus Ri is depicted for the TWRC with ideal/hypothetical AF relaying at SNR = 
20 dB when the sum rate R = Ri + R2 is, fixed to 0.5, 1, and 1.5, respectively. As can be seen 
from the figure, the bottleneck exponent E* (i?i, R2) is unimodal as a function of Ri for fixed 
R, and its maximum is at the mode of Ri determined by Theorem [2] for each value of R. We 
can also observe that the optimal/quasi-optimal rate allocation is of significant benefit to the 
bottleneck exponent. The optimal rate pairs (-Ri, -R2)opt sum rates R = 0.5, 1, and 1.5 are 

(0.176,0.324), (0.425,0.575), and (0.664,0.836), attaining the maximum E*{Ri,R2) equal to 
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1.0243, 0.5307, and 0.1995, respectively. For R > = 0.83, the quasi-optimal rate pairs at the 
sum rates R = 1 and 1.5 are (0.426,0.574) and (0.676,0.824), attaining E* {Ri, R2) equal to 
0.5286 and 0.1881, respectively. We can see again that the quasi-optimal rate pairs quite well 
approximate the optimal ones for R > R^ with a negligible loss in the bottleneck exponent. 

2) Optimal Power Allocation: Fig. [8] shows the bottleneck error exponent E* (i?i, R2) ver- 
sus R = Ri = R2 for the TWRC with ideal/hypothetical AF relaying under optimal and 
uniform power allocations at SNR = 20 dB. To determine E* (Ri, R2) in (|27l) . we first find 
the optimal power allocation p^^^ using Corollary [3] maximize i?"* (p^^^, p, Ri, R2) over p 
using the method given in [9, Section 2.2.4], and then successively perform the expectation 
of maXpg[o_i] (Popt7 P) -Ri) with respect to ctker by the Monte Carlo method. Compared 
with the uniform power allocation, we can see that the optimal power allocation significantly 
improves the bottleneck error exponent. 



VI. Conclusions 



In this paper, we have derived Gallager's random coding exponent to analyze the fundamental 
tradeoff between the communication reliability and transmission rate in AF two-way relay 
channels. The exponent has been expressed in terms of the generalized Fox iJ-function, from 
which the capacity and cutoff rate were also deduced for the link of each direction in the TWRC. 
Using the worst exponent decay between two links as the reliability measure for the information 
exchange, we put forth the concept of the bottleneck error exponent to effectively capture the 
tradeoff between the rate pair of the two links and the information-exchange reliability at such a 
rate pair for the design of two-way relay networks such that both links can communicate reliably. 
As its applications, we formulated the optimal rate and power allocation problems that maximize 
the bottleneck error exponent. Specifically, we presented the optimal rate allocation under a sum- 
rate constraint and its simple closed-form quasi-optimal solution that requires knowing only the 
capacity and cutoff rate of each link. The optimal power allocation under a total power constraint 
of the two terminals was further determined in the presence of perfect global CSI by solving 
the quasi-convex optimization problem. 
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Appendix 



A. Proof of Theorem [7] 

Let Vfe = pfcOife and = {PK + pm/Pk) OLXOL2I0L],. Then, 14 ~ £: (A^), W^^ ~ £^ (yUfc), and 



(30) 



Using the probability density function (PDF) of the Harmonic mean of the two exponential 
random variables [12] and the transformation py (y) = -^\Px (y/o) where Y = aX, we obtain 
the PDF of 7^^ as 



^7"^ (7) = ^Xkf^kie 
Tk 



\ Vk J 



+ — Afc + fik) V^'kl^kie "fc Ki , 7 > 

Vk \ Vk J 



(31) 



where (■) is the z/th order modified Bessel function of the second kind whose integral 
representation is given by [18, eq. (8.432.6)]. 
Using (|3TI) . we have 



Since it is obvious that i?o,fc (p) = for p = 0, we define 



p^ub (7) d-f } . 



(32) 



x(l + ax) e '"'K^ (cx) dx 



(33) 



to find Eo^k (p) in (|32l) for < p < 1. To evaluate the integral X (p), we first express (1 + ax) ^ 
and e^^Ki, (cx) in terms of the Fox i7-functions with the help of [19, eqs. (8.3.2.21), (8.4.2.5), 
and (8.4.23.5)] as follows: 



e^^'K^ (cx) 



r(p) 



ax 



cos VTt 



TT 



2,1 



(i-p,i) 

(0,1) 

(1/2,1) 
('^,1),(-'^,1) 



2cx 



(34) 
(35) 



where -ff^^" [■] is the Fox iJ-function [19, eq. (8.3.1.1)]. Then, substituting dH]) and ([35]) into 
(1331) . we have 
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COS [urr) 
0iT (p) Jo 



xe 



cos [uTr) 

v^r (p) 



-{b+c)a; tt1,1 
-"1,1 



/T , \-2 rjl,l,l,1.2 

(6 + c) ^i;(ia):o,(i; 



ax 



(1-P, 1) 
(0,1) 



rr2,l 
-"1.2 



2cx 



(1/2,1) 

('^, 1) , i-i^, 1) 



fix 



2) 



a 
b+c 

2c 
b+c 



(2,1) 

(l-p,l); (1/2,1) 



(36) 



(0,l);(z/,l),(-z/,l) 

where the last equality follows from [15, eq. (2.6.2)]. Finally, from (l3TI)-([33]) and ([36]) . we get 
(fT2l) and complete the proof. 

B. Proof of Corollary [7] 

It follows from Theorem [T] that 

(Ck) 



1 dEp^k (p) 

2 dp 

J p=0 

In (1 + 7) P ub (7)^7. 



(37) 



Similar to the derivation of Eq ^ (p), we first express In (1 + 7) in terms of the Fox //-function 
with the help of [19, eq. (8.4.6.5)] as 



ln(l + 7)=i72 



1,2 



7 



(1,1), (1,1) 
(1,1), (0,1) 



(38) 



Then, again using (|35l) and [15, eq. (2.6.2)], we evaluate (|37l) as (O and complete the proof. 

C. Proof of Theorem |2] 

Since the exponent E^^k (Rk) is a monotonically decreasing function in R^, it is obvious that 
for any rate pair ,^2) ^ ^ with the sum rate Ri + R2 = R> \Ro,i — Ro,2\, 

E:{R,,R2)>E:{R,,R2) (39) 

whenever R2) is such that i?r,i (-Ri) = -E'r,2 (-R2) and Ri + R2 = R. Therefore, the optimal 
solution (_Ri, -R2)opt of the problem Vi for R > |_Ro,i — -Ro,2| is uniquely given by 

i?2)opt e {(i?l, i?2) e ^ : (i?l) = ^r,2 (i?2) , i?l + i?2 = r} . (40) 

Although, clearly, the optimization problem (|20|) is mathematically challenging, it follows from 
(|40|) that the optimal solution {Ri,R2)opt — l-^o.i — -Ro,2| is the intersection point of the 

rate-pair curve ^ and straight line .if, and we can determine it graphically, as shown in Fig. HI 
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Let = {{Rl, R2) e ^ : Q < Rl < Rcr,i, < R2 < Rct,2} and be the largest sum rate at 
which the optimal solution (i?i,i?2)opt of the problem Vi belongs to the subregion ^i. When 
the rate is less than the critical rate, the optimal value of p is equal to 1 and the RCEE for the 
link LfcgT- of the TWRC can be written as 

Ei,k {Rk) = Eo^k (1) — 2-Rfc 

= 2 {Ro,k - Rk) . (41) 
Therefore, for the sum rate R < R^, the problem Vi can be rewritten as 



Pi 

which is equivalent to 



max miiifcgr {Ro,k - Rk) 
Ri,R2 

s.t. Ri + R2^R (42) 

< i?l < -Rcr,l, < i?2 < Rcr,2 



max min {Rq 1 — i?o 2 — R + Ri} 
Vi^{ ' ' (43) 

s.t. < i?i < R < i?cr,l + Rcr,2- 

Without loss of generality, we assume i?o,2 > Ro,i, and we can consider two different cases as 
follows: 

• When R < it!o,2 — Rq,i, we have it!o,2 — R + -Ri > -Ro,i — Ri and 

min {i?o,i - Rl, Ro,2 -R + Ri} 

— Ro,i — Rl 

< Ro,i. (44) 
Thus, in this case, the optimal rate pair is 

(i?i,i?2)opt = (0,R). (45) 

• When R > i?o,2 — Ro,i^ we need to consider two additional cases. 
If i?o,i -Ri> Ro,2 - R + i?i or i?i < i (R - Ro,2 + Ro,i), then 

min {i?o,i - Ru Ro,2 -R + Ri} 

= Rq,2 — R + Rl 

_ _ R — Rn 9 + Rn 1 
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Therefore, the optimal rate pair is given by 

/ D n X / + -^0,1 — Ro,i R ~ -Ro,i + Ro,i 

l^i,^2jopt- 1 2 ' 2 

If i?o,i -Ri< Ro,2 - R + i?i or i?i > I (R - i?o,2 + Ro,i), then 

min {i?o,i - Ri, Ro,2 - R + 
= -Ro,i — Ri 

^ D R — Ro.2 + -Ro,i 

- 2 

_ — R + i?o,2 + -Ro,i 

~ 2 

Therefore, the optimal rate pair is given by 

R + Ro,i — Ro,i R — Ro,i + -Ro,i 



(-Ri, R2 



'opt 



(47) 



(48) 



(49) 



2 ' 2 

Since -R2)opt should belong to Mi, we can find the decisive sum rate as (|22|) from 
the fact that 

R + Ro,i — Ro,i ^ D 

S -Kcr,l 

(50) 



R — -Ro,i + R^ 



0,1 



< -Rcr,2- 



From (|45I) . (|47I) . and (1491 ), we arrive at the desired result (|2i 



D. Proof of Theorem |3] 

For any t G IR+, the upper- level set of (i?^) that belongs to S is given by 

1 V'iPRttltt2 



In 



1 + 



+ • „/,2 



-p 



- 2pRk > t 



V^fcttfc + (pr + i'li'll'^l) Oiia2lak + 1 



> V 1 + PR"i"2/afc + \\Atpf 



1 +PRaia2/afc / . 



(51) 
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with 

A = diag {^/a^, \f(^ (52) 



A (1+P) 

Vk = 



exp 



P 



(53) 



where >ix. denotes the generalized inequality with respect to the second-order cone (SOC) /C 
[16] and is a standard basis vector with a one at the /cth element. It is clear that \J 
is a convex set since it can be represented as an SOC. Since the upper- level set U {E^^\.^£) is 
convex for every t E IR+, -E™^ (p, p, Rk) is, thus, quasi-concave^ 

We now show that E^^}: {p,p,Rk) is not concave by contradiction. Since the function In(-) 
is a monotonic function, we simply need to show that fk ("0) = ? , t'', o\ ; is 



not concave. We consider and such that = (Ck and tpf^ = 6(ek for < ( < V P and 
< 6 < 1. Clearly, and are feasible solutions of V2- For any A G [0, 1], we have 



-1 

A 



fk {M>a + (1 - A) ^fe) = (^ttfc + ^2 J = ^^^^ 

where (C) is clearly convex in C^. Due to convexity of {Q, the following inequality must 
hold 

Qk (ACa + (1 - A) Q < \gk (Ca) + (1 - A) gu (G) ■ (55) 

Now, by letting Ca = C/ (A + 5 (1 - A)) and Qb = 6C/ (A + 5 (1 - A)), we can rewrite ([55]) as 

h {\^, + (1 - A) ^,) < A/fc J + (1 - A) h M ■ (56) 

Thus, we have showed that there exist ipai'^b ^ ^+ A G [0,1] such that (l56l) holds. By 
contradiction, (V') is not a concave function on M^. Therefore, it follows that E^^^ (p, p, Rk) 
is also not concave. 

Since the nonnegative weighted minimum of quasi-concave functions is quasi-concave [16], 
-^r* {Pi P-iRiiR^) is also quasi-concave. Furthermore, V2 is a quasi-concave optimization prob- 
lem since the constraint set in V2 is convex in 



' Note that a concave function is also quasi-concave. 
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Fig. 2. Random coding error exponent for the link L^gr of the TWRC and OWRC with ideal/hypothetical AF relaying, 
ni = 0.5, ^2 = 2, and SNR = 20 dB. 
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Fig. 3. Capacity and achievable sum rate versus SNR for the Unk L^gr of the TWRC and OWRC with ideaI/h5rpothetical AF 
relaying. Q,\ = 0.5 and ^2 = 2. 
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Fig. 4. Cutoff rate versus SNR for the link L^gr of the TWRC and OWRC with ideal/hypothetical AF relaying. Qi = 0.5 
and Q2 = 2. 
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Fig. 5. Bottleneck error exponent E* {R\,R2) versus Ri for the TWRC with ideal/hypothetical AF relaying at R2 = 0, 0.2, 
0.5, 0.7, and 1.1 nats/s/Hz. fii = 0.5, ^2 = 2, and SNR = 20 dB. The values of minjiJi € E+ : Er,i (Ri) < £r,2 (-R2)} are 
equal to 0, 0.16, 0.36, 0.54, 0.92 nats/s/Hz for R2 = 0, 0.2, 0.5, 0.7, and 1.1 nats/s/Hz, respectively (indicated by the cross 
marks). 
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Fig. 6. Optimal rate pair (i?i,_R2)opt tliat maximizes tiie bottleneclc error exponent E*{Ri,R2) for tiie TWRC witli 
ideal/iiypothetical AF relaying at sum rates R = 0.148, 0.4, 0.8, 1.2, 1.6, 2.0, and 2.4 nats/s/Hz. fii = 0.5, ^2 = 2, and 
SNR = 20 dB. For R > = 0.83, the quasi-optimal rate pairs ai'e also plotted for R = 1.2, 1.6, 2.0, and 2.4 nats/s/Hz. 



NGO et at: TWO-WAY RELAY CHANNELS: ERROR EXPONENTS AND RESOURCE ALLOCATION 27 




Fig. 7. Bottleneck error exponent E* (_Ri,i?2) versus Ri for the TWRC with ideal/hypothetical AF relaying at sum rates 
R — 0.5, 1, and 1.5 nats/s/Hz. Q,i — 0.5, 02 ~ 2, and SNR — 20 dB. The optimal rate pair (Ri, R2)^p^ for each sum rate and 
the quasi-optimal rate pairs for R > = 0.83 are also plotted. 
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Fig. 8. Bottleneck error exponent E* {Ri, R2) versus -R for the TWRC with ideal/hypothetical AF relaying with optimal and 
uniform power allocations. Ri = R2 = R, fti = 0.5, ^2 = 2, and SNR = 20 dB. 
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Fig. 9. Graphical interpretation of the optimal rate pair {R\,R2), 



